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Abstract: Using a semiclassical approach to Gravitoelectromagnetic Inflation (GEMI), we study 
the origin and evolution of seminal inflaton and electromagnetic fields in the early inflationary 
universe from a 5D vacuum state. We use simultaneously the Lorentz and Feynman gauges. Our 
formalism is naturally not conformal invariant on the effective 4D de Sitter metric, which make 
possible the super adiabatic amplification of electric and magnetic field modes during the early 
inflationary epoch of the universe on cosmological scales. This is the first time that solutions for 
the electric field fluctuations are investigated in a systematic way as embeddings for inflationary 
models in 4D. An important and new result here obtained is that the spectrum of the electric field 
fluctuations depend with the scale, such that the spectral index increases quadratically as the scale 
decreases. 
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1. Introduction 

The origin of cosmological scales magnetic fields is one of the most important, fascinating and 
challenging problems in modern cosmology. Many scenarios have been proposed to explain them. 
Magnetic fields are known to be present on various scales of the universe Q. Primordial large-scale 
magnetic fields may be present and serve as seeds for the magnetic fields in galaxies and clusters. 

Until recently the most accepted idea for the formation of large-scale magnetic fields was the 
exponentiation of a seed field as suggested by Zeldovich and collaborators long time ago. This 
seed mechanism is known as galactic dynamo. However, recent observations have cast serious 
doubts on this possibility. There are many reasons to believe that this mechanism cannot be 
universal. This is why the mechanism responsible for the origin of large-scale magnetic fields is 
looked in the early universe, more precisely during inflation|^, which should be amplified through 
the dynamo mechanism after galaxy formation. In principle, one should be able to follow the 
evolution of magnetic fields from their creation as seed fields through to dynamo phase characteristic 
of galaxies. It is believed that magnetic fields can play an important role in the formation and 
evolution of galaxies and their clusters, but are probably not essential to our understanding of large- 
scale structure in the universe. However, an understanding of structure formation is paramount to 
the problem of galactic and extragalactic magnetic fields]^, ^. 

It is natural to look for the possibility of generating large-scales magnetic fields during inflation 
with strength according with observational data on cosmological scales: < 10~^ Gauss[^. However, 
the FRW universe is conformal fiat and the Maxwell theory is conformal invariant, so that magnetic 
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fields generated at inflation would come vanishingly small at the end of the inflationary epoch. The 
possibility to solve this problem relies in produce non-trivial magnetic fields in which conformal 
invariance to be broken. 

On the other hand, the five dimensional model is the simplest extension of General Relativity 
(GR), and is widely regarded as the low-energy limit of models with higher dimensions (such as 
lOD supersymmetry and IID supergravity). Modern versions of 5D GR abandon the cylinder and 
compactification conditions used in original Kaluza-Klein (KK) theories, which caused problems 
with the cosmological constant and the masses of particles, and consider a large extra dimension. 
In particular, the Induced Matter Theory (IMT) is based on the assumption that ordinary matter 
and physical fields that we can observe in our 4D universe can be geometrically induced from a 5D 
Ricci-flat metric with a space-like noncompact extra dimension on which we define a physical 5D 
apparent vacuum. The vacuum we shall consider is very restrictive in the sense that we shall not 
consider any kind of charges, matter or currents on the 5D spacetime. In a relativistic framework, it 
can be expressed by the 5D null geodesic equations, which are only valid for massless test particles 
in 5D. However, observers that move with frames If^ = = (described by a constant foliation 
on the extra dimension), can see the physics described by the effective 4D energy- momentum 
tensor embedded in the 5D apparent vacuum, which is geometrically described by a 5D Ricci-flat 
spacetime. From the mathematical point of view, the Campbell-Magaard theorem[^ serves as a 
ladder to go between manifolds whose dimensionality differs by one. This theorem, which is valid in 
any number of dimensions, implies that every solution of the 4D Einstein equations with arbitrary 
energy-momentum tensor can be embedded, at least locally, in a solution of the 5D Einstein field 
equations in vacuum. Because of this, matter, charge and currents may be 4D manifestations of 
the topology of space. 

Gravitoelectromagnetic Inflation (GEMI)0 was proposed recently with the aim to describe, in 
an unified manner, electromagnetic, gravitational and the inflaton fields in the early inflationary 
universe, from a 5D vacuum. It is known that conformal invariance must be broken to generate 
non-trivial magnetic fields. A very important fact is that in this formalism conformal invariance is 
naturally broken. Other conformal symmetry breaking mechanisms have been proposed so farjsj. 
However, most of these are developed in the Coulomb gauge. In order to simplify the equations of 
motion for A'^ , in this paper we use simultaneously the Lorentz and Feynman gauges, to calculate the 
electric and magnetic spectral indices for the spectrums of these fluctuations taking into account 
the induced currents. The main contribution of this paper is the study for the spectrum of the 
electric field fluctuations in a systematic way as embeddings for inflationary models in 4D. This 
topic has been ignored in the literature. 

The paper is organized as follows: in Sect. II we introduce the 5D vacuum of the fields on a 
generic 5D Ricci flat metric, to obtain the equations for the vector fields using simultaneously the 
generalized Lorentz and Feynman gauges. Also, we impose a semiclassical approach to the vector 
fields. In Sect. HI we study the particular case of a 5D Ricci flat space-time for an extended de 
Sitter expansion. In Sect. IV describe the dynamics of the vector fields on an effective 4D de Sitter 
space-time, when we make a static foliation on the noncompact extra dimension, which is considered 
as space-like: ip = ipQ. We develop the equations of motion for the fields using a particular Lorentz 
gauge on the effective 4D de Sitter space-time. After it, we describe the dynamics of the classical 
and quantum fields, to finally calculate the evolution and spectrums of the inflaton, electric and 
magnetic fields. The conclusions are developed in the Sect. V. Finally, we have included two 
appendixes where we have developed respectively the details of the calculations for the modes of 
the electric field fiuctuations, and the spectrum for these fluctuations. 
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2. Vector fields in 5D vacuum 



We begin considering a 5D manifold Ai described by a symmetric gab = gta^ 
This manifold M is mapped by coordinates {a;°} 

dS^ ~ gabdx'^dx^ . 



5D tensor metric. 



(2.1) 



From the geometrical point of view, to describe a relativistic 5D vacuum, we shall consider that gab 
is such that the Ricci tensor Rab — 0, and hence: Gab = 0. To describe the system we introduce 
the action on the manifold M 



S = / d^'Xy/^ 



16nG 4^"^^' 



(2.2) 



where i? is the 5D scalar curvature on the five-dimensional metric (2.1) and Q"''' = F"** — 
jg''''VfAf, where the 5D Faraday tensor is = W'^A" - V=A* = d^'A" - SM''. We shall consider 
that the fields A^ are minimally coupled to gravity and free of interactions, so that the second term 
in the action is purely kinetic. 



2.1 Einstein Equations in 5D 

If we minimize the action respect to the metric we will obtain Einstein Equations in 5D. In this 
paper we shall use a semiclassical approach where the Einstein equations are expressed by the 
homogeneous component of the fields. This slightly differs from the one used by |^ in the fact 
that we don't need to renormalize the stress tensor, but at the cost of assuming a semiclassical 
behavior of the fields that rules out the dependence with the wavenumber in the calculation of the 
semiclassical Einstein equations 



Gab 



-SttGT, 



(0) 



ab 



where xj^^'^ = {Tab{A'^)). Notice that we use a semiclassical expansion of the vector fields 



A" 



A" + 5 A" 



(2.3) 



(2.4) 



where the overbar symbolizes the 3D spatially homogeneous background field consistent with the 
fixed homogeneous metric and 5A'^ describes the fiuctuations with respect to A'^ . In this sense 
when we perform the expectation value of the stress tensor, adopting the ansatz {5Ac) = 0, only 
will appear zero order and the second order t'^ in perturbations terms. The last corresponds 
to a feedback term and is related to back-reaction effects, which do not will be consider in this 
paper. The stress tensor is defined by the fields lagrangian being symmetric by definition 



Tb, 



2 f d 



d 



d 



dg 



be 



(2.5) 



The appearance of variations with respect to derivatives of the metric is because we are dealing 
with vector fields whose covariant derivative operators involve Christoffel symbols (i.e. ordinary 
derivatives of the metric). In our case the stress tensor reduces to 



Tbc — F'^bFce + -^gbcFdeF'^ 



- A 2A^ 



A, 



56c 



A' 



{b-c) 



(2^(fc5c)/.,/ 
3 



9hf., (b 



K) ) 



^Tl,TlfA''\Af 



gbcjAfA^,] 



(2.6) 



(2.7) 



where 7^ 



2A 

5 • 



^In our conventions latin indices "a,b,c,..,h" run from to 4, greek indices run from to 3 and latin indices 
"i,j,k,..." run from 1 to 3. 
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2.2 5D dynamics of the fields 

The Euler-Lagrange equations give us the dynamics for At 

WfWfA'>^R^j:Af-{l-X)W''WfAf = 0. (2.8) 

In particular, the choice A = 1 is known as Feynman gauge, somehow equivalent to a Lorentz 
gauge '^fA^ = 0. In this paper we shall choose simultaneously both conditions. The first one 
assures that the balance of each component of any external current to be null, and the second one 
is more restrictive, because assures that each direction (insider or outsider) of each component of 
the current to be zero. 

It is easy to show that the 5-divergence of the field equation of motions satisfy the same equation 
as in a Minkowski space, but changing ordinary partial derivatives by the covariant derivative 

V^V^V/A^) =0. (2.9) 

Hence, the Lorentz gauge is satisfied for appropriate initial conditions of Va^" = 0. With such a 
choice the field lagrangian density Cf = —jQ^ is 

C'f = -iv.AbVM" = -iv^A.V'A'' - iv4A.V4A'' - ^V^A^WA* - ^V4A4V^A\ (2.10) 

where — Q'^''Qab- For 4D observers living in a hypersurface where the fifth component of the 
vector field is normal to it, this extra dimensional field will manifest separately, like an effective 4D 
vector field A'^ and a 4D scalar field A*. In this sense we can identify kinetic terms for both, scalar 
and vector fields, and the derivatives with respect to the extra dimension may be interpreted as 
potential (or dynamical sources) terms joined with massive terms for each of them. 

The stress tensor in this gauge is 

Tab^ ~VaAeVbA''-VeAaV''Ab~2g,^,A,)TlfV'Af + ^gabVeAfV'Af ~ (2.11) 
2^ [V(aA)^-^' + y^A^bAa) - A(,Vb)A^] - [ViaAb)Af + V^A^hAa) - A(,Vb)A/] . 

9 '■' 

3. Special case: 5D generalization of a de Sitter spacetime 

Because we are interested to study a cosmological scenario of inflation from the context of the 
theory of Space-Time-Matter, we shall consider the 5D Riemann-flat metric |pO| 

where is a time-like dimension related to the number of e-folds, dr^ — dx^Sijdx^ is the Euclidean 
line element in cartesian coordinates and ip is the space-like extra dimension. This metric satisfies 
the vacuum condition = 0. 

For this 5D metric the field equations, after taking Lorentz gauge: WaA'^ ~ OnA'^ + 3A'^ + 
d^A'^ + Aij-^A^ + d,A' = 0, are 



r 52 
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ip dip 
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A^ - 2d^ {a° + — 
12 



^4 = 0. 



ip_ 
= 0, 



(3.2) 
(3.3) 
(3.4) 
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Notice that the (3.4) is decoupled after applying the Lorentz gauge. However we see that it is not 
sufficient to decouple all the field equations. This is because the non zero connections of the metric 
(|3.l|) act in a non trivial manner in the vector fields derivatives. There are 14 non zero Christoffel 
symbols 



-^4-'. , r^o = i. r°=e2^, r^,^^, = -^e^^. (3.5) 

Therefore, in this Riemann-flat spacetime we obtain the D 'Alambcrtian of the field 



VfV^A'' = 0, 



(3.6) 



but, expressed in terms of the ordinary derivatives and the Christoffel symbols we notice the coupling 
terms 

gf'^ {dfd^A" + 2VlfdhA^ + fA^ - T%d,A' - T),j:l^A'' + tIjTI^A''} = 0. (3.7) 

Notice that in a 5D Minkowskian metric: dS^ — dt^ — dr^ — dip'^, the connections vanish and the 
field equations remain decoupled after the gauge choice. 

3.1 Dynamics of the 3D spatially isotropic background fields 

We shall combine the field equations of motion for the classical homogeneous fields with the Einstein 
Equations, the first ones reduce to 



2 






2 






d 
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dN 


9^2 



92 6 9 

9?/'2 d'\\) 



+ 



6 d 



12 



AO 

=0, 

A^ =0. 



2 d 



A^ = Q, 



(3.8) 
(3.9) 
(3.10) 



■0 dtp ijP 

Notice that the equation for A^ is the unique coupled. Furthermore, once obtained A**, we can 



describe the dynamics of A^ in (3.8), where A"^ appears as a source 



4. Effective 4D dynamics of the fields 



Now we consider a static foliation on the 5D metric (3.1). The resulting 4D hypersurface after 
making = V^o describes a de Sitter spacetime. From the relativistic point of view an observer 
moving with the penta velocity — 0, will be moving on a spacetime that describes a de Sitter 
expansion which has a scalar curvature R = Vl/^^ = 12 H^, such that the Hubble parameter is 
defined by the foliation Hq = ip^ 



. Hence, if we consider the coordinate transformations on 

2 



we then arrive to the Ponce Leon metric||TT|: ^5*2 = (^^j [df - e2*/^"di?2] _ dip"^. If we fohate 
-0 = -00, we get the effective 4D metric 



dS-" 



ds' 



dt^ 



e^"°^dR^, 



(4.2) 



which describes a 3D spatially flat, isotropic and homogeneous de Sitter expanding universe with a 
constant Hubble parameter 7?o. 



The dynamics of the fields being given by the equations (3.2), (3.3) and (3.4), evaluated on 



the foliation ip = %pQ — I/Hq, with the transformations (4J_). In the following subsections we shall 
study separately the dynamics of the classical 3D spatially isotropic fields: A'^{t, ipo) and ^^{t, ipo), 
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and the fluctuations of these fields: SA'^{t,R,iljo) and 6A'^{t, Rjipo). Notice that now R = R{X^). 
To describe the dynamics of the fields we shall impose the effective 4D Lorentz gauge: V^^A'^ = 0. 
It implies that the 5D Lorentz gauge with the transformations (O) and evaluated on the foliation 
must now be 



0, 



(4.3) 



where '^^^V^ denotes the covariant derivative on the effective 4D metric (4.2). Hence, in order 
to the effective 4D Lorentz gauge to be fulfilled, we shall require 



{d^A^ + 4i;-^A^ 



= 0. 



(4.4) 



4.1 4D classical field dynamics 



In order to solve the equations (3.8), (3.9) and (3.10) on an effective 4D de Sitter metric, we must 
evaluate these equations on the particular foliation = ipo = Hq^, r = Ri/jq and N — Hot. 
We shall identify the effective scalar A'^ with the infiaton field: A'^{t, R,iPq) = (p{t,R,ipQ) 
and we shall denote (j){t,tpQ) ^ 't'li^) 4'2{'>ij)\j^^Hot ^=ipo=H^^^ spatially isotropic and 

homogeneous background field. In the same way we state for the homogeneous component of the 



vector field the separation A^{t,ipo) ^ S{{N)S2{ip) 
index j to label the functions Si{t) and S2{ipo)- Hence, we obtain 



^ , in the next we shall drop the 



Am? 



where we have considered the condition (|4.4|), such that 



dP_ i_d_ 



m'^ (l){t,'iJjo), 



(4.5) 



(4.6) 



4'o 



where (j) plays the role of the background infiaton field. Furthermore, the general solution of eq. 



( |3.9D on the effective 4D metric (|4J), is 



A^it,i,o)-Sit)=e- 



ci e ° + C2 e ° j , 



I - 



4i^ 
~2b 



where 



6_d_ 

ip dip 



A' 



(4.7) 



(4.8) 



A similar treatment can be done for after making use of the condition (4.4), the transforma- 
tions (4J_) and the foliation ■0 = ■0o = f/-ffo- However, the difference with the other background 
components of the field observed in eq. (3.8) is that A^ = (j)(t, ipo) acts as a source of A'^{t,ipo)- 

As a particular choice we shall consider a 4D inflationary universe, where the background flelds 
are A^ = (0,0, 0,0, i^), in agreement with a global (de Sitter) accelerated expansion which is 3D 
spatially isotropic, flat and homogeneous. ^. In this case, the relevant components of the classical 



■^One could consider, for instance, the case when the background field is = (<&, A^, 0, 0, O) , that defines an 
effective homogeneous component of the electric field. However, we would obtain an anisotropic component of the 
stress tensor Tio, which is not compatible with our background, spatially flat, homogeneous and isotropic (de 
Sitter) metric. In general this implies that for the background fields to satisfy Einstein equations, the components 
Ao\ Ai; A2\ As are highly restricted. In particular we have the following cases to choose:(i)A' = 0, A" = i/iq) 
and A** = (f>{t, i/iq), {ii)A'' = A* = and A' = Aq constants. In what follows we shall analyze a particular choice of 
the first case (with A" = 0), because the other isn't very interesting in the physical sense. 
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Energy momentum tensor, are 



1 -2 



{Tp)\a=tl3 = 0, 



5 72 












^ ^2, 




- —4>4>' 









(4.9) 

(4.10) 
(4.11) 



where the prime denotes the partial derivative with respect to and dots denote partial derivatives 
with respect to the time, which in our case are zero: 4> — 0. Furthermore, from eq. (4.9) we can 



make the following identification for the background scalar potential: 



V\<h] = 



5 



2 



(4.12) 



In our model, the hypersurface tJj = t/jq defines a de Sitter expansion of the universe with a Hubble 
parameter Hq = i^o^- The equation of state for this case is p = —p= —3/ (SttGiPq). Then, it is 
easy to see that the only compatible background solution for the field evaluated on the hypersurface 
is the typical de Sitter solution for a background scalar field: (i>{t = N/Hq,^o) = (f>o. This means 
that 



A particular solution of (4.4) is 



-4:HQ(j)Q. 



it>=->Po = l/Ho 



From eqs. ( 4.12 ), ( 1.13 ) and the second in ( 1.14 ), we obtain 



(4.13) 



(4.14) 



(4.15) 



(4.16) 

It is easy to see by inspection in ( |3.8[) that <t>{t, tp) is a constant of N. In other words, the unique 



4,=,f,„ = l/Ha 



40 7rG" 



such that replacing (4.15) in the second equation of (1.14), we obtain 

' 57rG' 




origin of the effective 4D potential energy density (4.13) related to the background infiaton field is 
the ■0-dependence of (f>{N, tp). 

4.2 4D Field fluctuations 



Here we consider equations (3.2), (3.3) and (3.4) to search for possible electromagnetic fields gener- 
ated through this model. In Sect. (4.1) we've seen that the Einstein equations for the background 
fields exclude any possibility of spatially homogeneous electromagnetic fields. 



The equation for the efi^ective scalar 6A*{t, R, t/jo) on the effective hypersurface (4^) is decoupled 
from the dynamics of the 4- vector. In contrast, the equations for SA^{t, R,ipo) and SA''{t, Rjipo) 
remain coupled. By the use of our 5D Lorentz gauge evaluated on the foliation = ^po = Hq^: 



i>o=H„ 



= 0, we can express the inhomogeneous term for 5A^ as only a function of SA'^. 



The solution will involve both, homogeneous and inhomogeneous parts. Once obtained 5A^ and 
SA^, we can finally search solutions for the components 5A^ . These total solutions are necessary 
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to deduce the effective electric fields. In contrast, as we previously said, the equation of motion for 
pure magnetic fields may be obtained by just applying the curl in the 3-space to equation (3^). The 
last term in (3^) vanishes because is a 3-gradient, and so magnetic fields equations are decoupled. 
To quantize the field fluctuations on the effective 4D de Sitter spacetime (4.2), we shall consider 
the equations (3.2), (3.3) and (3.4), with condition (4.4), the transformations (4.1) and the foliation 
ip = ijjQ = I /Hq . The equal time canonical relations are 



M,(i,i?,Vo),ff(i,-R',V'o) 



(4.17) 



tlJo = l/Ho 

where g*-' are the space-like components of the tensor metric in (|4.2| ) and 6^^^{R — R') is the 3D 
Dirac's function. Furthermore, the canonical momentum is given by the electric field IP = = 
"V^ — V°A"'. The equations (3.2), (3.3) and (3.4) with the transformations (4.1) can be evaluated 
on the foliation Tp — ipo — ^ /Hq to give the dynamics on the effective 4D spacetime ( |4.2| ) . If we take 
into account the conditions (4.4), the effective 4D dynamics of the fluctuations describe an effective 
4D Lorentz gauge, so that 



5Ho- 



dt 
dSA^ 



+ ly^H^SA" = -2i?o'^, 



dt 



Hle-^"^^'dl6A^ 



H^SA^ = 2Hld^ {5A° + Hat 
56 = 0. 



(4.18) 
(4.19) 
(4.20) 



describe the 4D dynamics of the fluctuations. A very important fact is that the electromagnetic 
field fiuctuations 5A^^ obey a Proca equation with sources. 
The expansion of the free field in temporal modes is 

3 



^ ' A=l 



-iK-R 



t 

{k,\y 



JKR 



(4.21) 

The equation of motion for the temporal modes W{K,t,^Q) of the free contravariant fluctuations 
SAi" is ^ 

^ + 5i/o^ + [K^e-^"«' + i^^H^] I = 0, (4.22) 

where K = Hq fc (fc is a dimensionless wavenumber). Furthermore, £^{k, A) are the polarizations , 
such that in the Lorentz gauge the following expression holds; 

3 / 



(4.23) 

A=l \ / 

where we have introduced the effective mass rn^jj: = Hq{v^ — ^) of the redefined temporal modes 
Uxit) = e^"°''^'^W{K,t,'ipo), that obey the harmonic equation Uk + uj\{t)UK = 0. The time 
dependent frequency is defined by the relation K^K^^ — ni\ff. 

^l{t) = K// + {e-"^'Kf] ■ (4.24) 

1 /2 

Modes with uij^ > are stable, but those with uj^ < [i.e., with k < (25/4- i/'^) e"°^], are 
unstable. In the small wavelength limit these behave like plane waves in Minkowski space. Fur- 
thermore, the annihilation and creation operators a(^K.\) and comply with the commutation 
relations 



t 



(4.25) 



^parenthesis denotes that sum do no run over these indices. 
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25 



x{t) = — e 

^0 



Hot 



The solutions for the temporal modes is 
W{K,tM) = e-5««*/2 {ciHS^) [x{t)] + C2HS^) , 

(4.26) 

fl 2') 

where ■H^ '^[a;(t)] are the first and second kind Hankel functions respectively. We can also obtain 
the temporal modes for the covariant 5A^ which are related to the contravariant ones: TKif) = 
g2ffot W{K, t, ipo). The commutation relations (4.17) yield the following conditions over these modes 



-5) {TKn-V<fK) = ^&)e 



-Hat 



'eff 



'eff 



From these relation we can deduce the following apparently independent equations 



Q-Hot 



(4.27) 

(4.28) 
(4.29) 



which are only valid on (short) wavelength modes for which ojj^ > 0. Equations ( 4.28 ) and ( 4.29 ) 
give us the normalization conditions for the modes of (5A^t. On the other hand, these modes are 
unstable on cosmological scales: tuj^ < 0, and the expression (4.28) tends to zero. To apply these 
conditions wc take the very small wavelength limit for the Hankel Functions x{t) ^ |cr^ — i|. These 
means that K/Hqc ^ ''^eff^ WK{t) — Ke^^°*'. In this limit the conditions (4.28) and 

( 4.29|) become dependent one of the another, since 



-Hat 



K\UV 



2wK{t) 



uv 



1 

2K' 



Letting us choose ci = (Bunch-Davies vacuum), the solution for the modes is 



(4.30) 



4.2.1 4D electromagnetic fluctuations 

The electric field for a observer in 4D is defined by its 4- velocity = F^xu'^. If we choose the 



particular co-moving frame 

E,, 



(ffoV'o)-\0 



we obtain 



0, 



dX 



^2Hatd^^, 

dt 



2Ho e^"°*6A 



(4.31) 



The magnetic fields are defined by = \evXafiU^F°'^ , where e^Xap = \jY'^''g\A,jXai3 is the totally 
antisymmetric Levi-Civita tensor and AyXap is a totally antisymmetric symbol with y^oi23 = ^1- 
Then for a co-moving observer we will have a magnetic field. 



Bn = 0, 



Tj^kl 



AjQki F 



Prom the last expression we can arrive to another that will be useful to obtain an equation of 
motion for the magnetic fields, we first define the Levi-Civita symbol in the 3-flat space using the 
co-moving frame: ejki = Ajoki (we note that 6123 — 1). Hence 



B, 



(4) 



9 9' 



^^'u^'e.udk'A'. 



(4.32) 
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For our particular case we obtain 



5^^ ejki dk' 



A' 



(4.33) 



The differential operator between square brackets commutes with the one applied to A' in the 
equation ( |3.3D , so that in the equation of motion for Bj — e^^^^Bj there will be no sources. We 
can express the field in Fourier components of the SA' field 



Here Vj{K, t, -ipo) — —iKj W{K, t, ipo) are the temporal modes with their complex conjugate V*{K, t, ipo) 
iKjW*(K,t,'4!Q). We perform the vacuum expectation value of the B-ficlds quadratic amplitude, 
defined by the invariant product (_B^) = (0|i?"i?Q|0). For comoving observers B'^ = and so we 
have = B^ Bj = e-^""* J^j -S/ J2j ■ Then 



(2^ 



(4.35) 



We will cut the above integral up to wavelengths that remain well outside the horizon wavenumber 
kn — cre^"*. In this limit we use the asymptotic limit of the Hankel functions for the long wavelength 
limit fce^^"* <C ^/a + 1. The power spectra is then 



2^"r2(g)g4 

47r3 



(4.36) 



if we ask for an almost scale invariant spectrum, then = | + ry, r; 



quadratic amplitude is then 



(B' 



50-^ 

y 



and < 1. The 



(4.37) 



where ^ 1 is a control parameter, such that we stay with super Hubble wavelenghts: k < Okn- 

Using only the homogeneous solutions of the equations ( 4.18| ) and ( 4.19D we can deduce their 
contribution for electric fields on the infrared (IR) sector, we obtain for comoving observers {E'^) m = 
{El +El + El)iR, where 



{E\)ir ^ -H, 



-AHot 



{E 



BIIR 



{Ec)iR 



° - f ) 



^-2H„t 



ek 







" dk 
2^ 



dk 
2^ 



Be 



2Hot , ^O"- 



'eff 



:{-^Hok,){Tkn-n% 



(4.38) 
(4.39) 
(4.40) 



The corresponding power spectrums are 



VEsik) 



22'^r2(a)ff4 



87r3 



{a^+a+ 1/4) ( 3e(-l+2-)«otfc3-2. ^ ^-2 ^2a Hot ,^5-2a 



VEcik) = 



(4.41) 

(4.42) 
(4.43) 
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The last goes to zero in cosmological scales since it is proportional to the wronskian (4.28). If we 



choose a ■ 



+ r], T] = -i^ and \v^\ < 1, we get 



{E\)iR 



BllR 



36*- 



4 
25 



{El 



C/IR 



0, 



(4.44) 

(4.45) 
(4.46) 



on cosmological scales. Notice that {E"^) is not scale invariant for a scale invariant magnetic field. 
Then we can say that on very large scales the amplitude of electromagnetic fields are 



3\/5 
27n7' 



1 /9 35/2 

kE')Tn - t ^0^""*^"^ 



(4.47) 



which are related to comoving observers. During inflation, the strength of the magnetic field in a 
physical frame is 



1/2 



-2Hot 



ot ^^2^ 



2\l/2 



IR 



(4.48) 



1 /2 

where {B^) ^ is given by the first equation in ( 4.47 ). At the end of inflation (i.e., for t — te), the 
size of the horizon was close to 3.6 x 10^^ cm. It has suffered an exponential growth ~ 4.4 x 10^^ (we 
suppose that the number of e- folds is = 63), from its initial value at Planckian scales. Hence, 

/ 2\ 

we can make an estimation for the strength magnetic fields at the end of inflation / 
on cosmological scales 



IR 



1/2 



IR 



3%/5 



X 10" 



(4.49) 



10^44 Ml ~ 10^6 Gauss, where 



0.223 X 10^" Gauss (I 



If we take Hq = 10^^ Mp, it holds 
Gauss ~ 0.6476 x 10^^^ GeV^)''. However, must be noted that this value is very sensitive to the 
number of e-folds suffered during inflation. 

On the other hand, the present day size of the universe is of the order of 10^® cm. (for 

/ \l/2 

t = to = 1-26 X 10^1 " ). We shall suppose that, after inflation ( B'^f^^^ j decreases adiabatically 

as a^^: {a{to) / a{te)) — 10^^^, so that the present day value for residual magnetic fields should 
be of the order of ~ 10^^^ Gauss. Of course, in this estimation is omitted any possible mechanism 
for the amplification of these magnetic fields [0, which could be taken into account. 

4.2.2 4D inflaton fluctuations 

For the fluctuations of the inflaton fleld we can make a similar treatment. The Fourier expansion is 



(t, R, Vo) 



(27r)3 



a(^)0(if , t, Vo) e'^-« + {K, t, Vo) e 



-iKR 



(4.50) 



such that the annihilation and creation operators cx(k.\) ^^id q;|^ ^-^ , comply with the commutation 
relations 



= (27r)^ S^^\K-K'). 



(4.51) 



In all the paper we consider natural units: h - 



- 11 - 



The solutions for the modes (j>{K,t,iJjo), are 

^{K,t,^lJo) = e-'"'''/^ {ciJ^[xit)]+C2Y^[xit)]}, ^=^^^m^. (4.52) 

The nearly invariant spectrum of the scalar perturbations is obtained for small values of the effec- 
tive mass: |m^| <C 1. After normalization of the modes, we obtain the standard result (see, for 
instancepH), on cosmological scales 



with amplitude 



«fc,(,*o) = Jj^«™r"°'l, ^=^--m' (4,53) 

which is divergent for an exactly scale invariant power spectrum corresponding to a null value of 
the inflaton field mass m. 

4.3 Effective 4D electromagnetic fluctuations with sources included 

In this section we shall find inhomogeneous solutions for the Fourier components of the fields, 
we have noted previously that magnetic fields are only generated through homogeneous solutions. 
Instead, electric fields are affected by the coupled dynamics of the equations of the model. This 
couplings come from the 5D background, because some connections in the 5D metric are not null 
[see (3.5)]. The equations (including sources) (4.18) and (4.19) may be written as 

1^ + 5i7o^ + e-'^«* + J/'i?o'| X^'^T^', (4.55) 
with an inhomogeneous solution 

X'^{t,k,^Po) ^""^ I dTF^'{T)ei"°^ {J^[x{T)]J^,[x{t)] ~ J,[x{t)]J^„[x{T)]}. (4.56) 

2Hq sm(o'7r) J 

where J^^{t) are different sources for each of the equations. Using the identities for the Bessel 
functions and their derivatives we arrive to the following expression for Fourier transform of the 
source term of ( 4.18 ): 

F\t) = -y^e-i^«* {(3/2 - nWi^Mt)] + x{t)Hf_Mt)]} ■ (4.57) 

Notice that the sources were omitted in a previous treatment |l^ . However, such that sources 
should be important, mainly for electromagnetic fields. 

Once known the solutions [see appendix(^)], we can define the Fourier components of the 
electric field: 

El^^{K,t) ^ ^iK^ X^{K,t) - Xl{K,t) - Xi{K,t) - Xl{K,t). (4.58) 

Here, the sufiix {X) means that we are dealing with the electric field calculated only with the 
inhomogeneous contribution of de modes in A'^: X'^{k,t,'4'o)- 

The amplitude of these fields on cosmological scales [i.e., the infrared (IR) sector], is given by 
the expression 

^w)^ i, j2^^T.Hx)iK.t)EU;{K,t), (4.59) 
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which has a power spectrum 

^^W = ^El^(x)P- (4-60) 
using the solutions we may write the power spectrum in the aproximate form [see appendix (^] 

oo 

Pk{t) ^k^J2 aq{ke-"°'Y°+'i (4.61) 

the dominant contribution for the electric field comes from the smaller spectral power with q ~ 0, 
setting like previously /Lt = 3/2 + e and (T = 5/2 + rywe obtain 

£;2 \ ^ e^^»* ^" , , (4.62) 

The first correction comes from 02 coefficient, since there is no ai. When we consider this term, we 
have no longer scale independence of the spectral index, and 

^n..-(^. (4.63) 
ao in k 



We shall only write an approximated expression for 02. As it is shown in ( |B.9[ ), X| dominates, so 
after considering only this contribution, one obtains 

77492 ,2 

where — ~ — | . The first correction due to the inhomogeneous contribution of the modes of A'^ to 

ao 3 <= 

the electric field amplitude is 

^fi)>-(^fo))(l^i^^), (4.65) 



where we remember that 9 = k/kn <C 1, kn = cre^^*" being the wavenumber related to the Hubble 

1/2 / vi/-^ 

radius in a comoving frame. In a physical frame we obtain (i?^) ^ a? (E'^hysf and the energy 
density: p ^ a^Pphys- Notice that the energy density related to the electric fields at the end of 
inflation {Hq ~ 10~^ Mp), is very small with respect to the background inflaton energy density: 



El A iE'^ 



I phys \ / phys 



10-2 ( ^ 



P V 

so that back-reaction effects due to the electric fields are really negligible during inflation. 

The figure (|^) shows Snk as a function of 0. Notice that Suk decreases almost quadratically as 
the wavelength decreases. When the horizon entry (i.e., after infiation when 9 = 9^, = k^,/kH, = 1), 
the value of Suk is close to Suk, — —0.035. However, during infiation the cosmological scales 
wavenumbers are k/ku < 10~^ that corresponds with Suk > — 10~^. Notice that we have taken 
into account the value kjj^ — ae^^ as the wavenumber related to the horizon wavelength when, 
after infiation, the horizon entry. 

5. Final Comments 

We have shown how primordial electromagnetic fields and inflaton fluctuations can be generated 
jointly during inflation using a semiclassical approach to GEMI. We have used simultaneously the 
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Lorentz and the Feynman gauges. The first one assures that the balance of each component of 
any external current to be null, and the second one is more restrictive, because assures that each 
direction (insider or outsider) of each component of the current to be zero. This is done with the 
aim to assure a 5D vacuum on the 5D Ricci flat metric (^^). In correspondence with this concept 
of vacuum, we have defined a 5D totally kinetic Lagrangian density £/ = —jQ^, which is totally 
absent of any kind of interactions. 

One of the important facts is that our formalism is naturally not conformal invariant on the 



effective 4D metric (4.2), which make possible the super adiabatic amplification of the modes of the 
electromagnetic fields during infiation in a comoving frame on cosmological (super Hubble) scales. 

In this paper we have analyzed the simplest nontrivial configuration field: A'' = [O, 0, 0, 0, 0(i, V'o) 
For this configuration of the background fields, the background infiaton field must be a constant on 
the metric (4.2) to satisfy the Einstein background equations in a de Sitter expansion: (f){t, ipo) = 4'o- 



Then, in the model here developed, the expansion of the universe is driven by the background 
infiaton field 0o and background electromagnetic fields are excluded to preserve global isotropy. 
Notice that back reaction eff'ects are not included, because the EM field does not contribute to the 
background expansion of the universe [Q, but however comes into play an important role at the 
perturbative level as vectorial metric fluctuations which are the geometrical reaction to the vector 
physical fields [p^. 

To describe the efi^ective 4D dynamics of the fields, we impose the effective 4D Lorentz gauge 



'^^^V^^'^ = 0, given simultaneously by conditions ([f.5|) and (4.4). Therefore, the origin of the 



generation of the seed of electromagnetic fields and the infiaton field fiuctuations during inflation 



can be jointly studied. The dynamics of 5A^ on the effective 4D metric (4^) obey a Proca equation 
with sources where the effective mass of the electromagnetic field fiuctuations is induced by the 
foliation ip = ipa = I/Hq- From the point of view of a relativistic observer this foliation imply that 
the component of the penta-velocity — Js ~ ^■ 

We have obtained that for small values of v a nearly scale-invariant long wavelengths power 

1/2 

spectrum for (-B'^) , which grows as a during infiation on a comoving frame. However, on a 
physical frame it suffer a super adiabatical evolution, so that at the end of inflation is of the order of 
10^^ Gauss. After inflation we have supposed that the fleld evolves adiabatically as a~^, to estimate 
the present day values on cosmological scales (for a physical frame): (-Bohus)^^^ — 10~^^ 

^ " Now 

Gaussjlj. On the other hand, the dominant terms in the amplitude of (E'^) grows as on a 
comoving frame, and has a scale dependent power spectrum with a spectral index ~ 3 + Suk- 
This is the main result of this paper. This scale dependence is described by Suk, which decreases 
quadratically as the scale decreases. In the limit case where fc — (very large scales), one finds 
that Snk=o — > 0. Finally, in what respect to the infiaton field fiuctuations {ScjP'), we obtain that 
they are nearly scale invariant on cosmological scales, and the amplitude is freezed in agreement 
with the predictions of standard 4D infiation. 
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A. The modes of the electric field 



In order to solve the integrate in ( [I.5q ) we express all the Hankel functions in terms of the first 
kind Bessel functions Ja[a;(i)], and J-a[x{t)\ 

hW(:,)= '^-"^-'"'"""'^"^'^^ (A.l) 
I sm(tT7r) 
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J—a(x) 6 Ja{x) 



—i sin((T7r) 

and then we expand them in their series representation 

(—1)'" 

m— 

such that the product identity isp5[| 

(-1)™ r{l + 2m + a + P) 



2m-\-a 



m=0 



m\ r{l + m + a)T{l + m + /3)r{l + m + a + 



2 



The terms included in expression (4.56) are of the form, 



dte"°\j^Jp)[x{t)]^—Y,' 

^ m,n 



mm 



x{t) 



2(m+n)-l+Q+^+7 



^ 

!n! '""■'^''^ I 2 



2(m+n)+Q+;3+7 



(A.2) 
(A.3) 

(A.4) 

(A.5) 
(A.6) 



where the coefficients C™^ ^ and 2?™^ ^ are defined by the following relations of the Gamma func- 
tions 



r(l + 2m + a + l3) 



(1 - 2m - a - ^)r(l + m + a + /3)T{1 +m + a)r{l + m + /3)r(l + n + j) 

r(l + 2m + a + /?) 
(-2m - a - /3)r(l + to + a + /3)r(l + m + a)r(l +m + /3)r(l + n + 7) ' 



, (A.7) 
(A.8) 



After some algebra, one arrives to the inhomogeneous solution for the modes of the electromagnetic 
field A° 



^3 /2\i^(-i)™+« /x(t)y(™+"' ' 



sin((T7r) y Ho \k 



E 



Ecrar 



2 



(A.9) 



where x{t) = ke The sum over s goes through three different powers. The coefficients £™ 

depend on the parameters fi, a and sum on indices m, n in the following way 



mn 



sin(/i7r) 



^ ^ ,, cmn ^ ^2 ( nmn nmn \ 



2 sm(/^7r) 



sin(/i7r) 2 



(A.io) 

(A.ll) 
(A.12) 



The inhomogeneous solution, {k,t,ipQ), of , has essentially three contribution terms. For 
simplicity, lets split the sources in the following way: 

J'i^2H^{zk^)W^°\k,t,i^o), Ti^2{iy)Hl<j,{k,t,i>o), J'^^2H^itk^)X^°\k,t,i^o)- 

Hence, the final solution of ( 4.56 ), written as X^ = Xf + X^ + X^, after using k^ — ke^ [e^ being 
an unitary vector), is given by the expressions 



X^(A,i,Vo) 



3in^(cr7r) V Hq \k 



7r3 ^2\^^ (-1)" 
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E 

m—l 

oo 

E 



r(i-a)r(i + a) 



r(l + n-cr) r(l + n + cr) 



a,^ (2) 



X \ f +"■+2™' 



(2) (^-^^^ 

, (A.15) 



with 



sin(cr7r) \J Hq \k 



n = 3/2 + ^, 



min 



^mn 



sin(/i7r) 



/ nvan n \ 

V'-'/j,(T,-cr '--/j,-<T,o-j 



1 



r, = 3/2 - 5™" = (C-,_. - C-_.,.) . 

Finally, the contribution of the inhomogeneous source is 



(A.16) 
(A.17) 



r3 /2 



47r 



sin((T7r) y i^o V^/ sin(o'7r) 



E 



nm!n!/i!r(l + « + c7)r(l + ;i-(7) V 2 



2{l+m+n+h) 



m + n) - |] - cr2 



2 



B. Calculation of the spectrum for the electric field fiuctuations 



(A.18) 



It is important to notice that X{ has appreciable differences with the other terms in ( 4.58| ); the 
one has a preceding factor fc^^/^, while the others {iK^X^, and X3) are proportional to fc^^/^. 
Furthermore, X{ has terms with the linear factor Hi^t while the others doesn't. 

Taking into account the last observation we arrange the power spectrum as follows 
00 00 00 

Vk{t) = fc' E aq{ke-"°'f°+'i + fc ^ bq{t){ke-"°')^°+'i + ^ c,(t)(fce-^°*)''«+«. (B.l) 

9=0 q=0 9=0 

The coefficients aq are all calculable from quadratics and cross products of: iK^ X'^ , and ATg. 
The coefficients 6,(t) = h^q^ + b^qh are linear in time and come from products of X{ with the 
others. Finally, the coefficients Cq{t) = c^q '' + c^qH + Cg'^'i^ are quadratic in time and are found 
from lAT^Jp. The lowest powers from which each term in the series begin are: /3o = 3 — 2/i, 
7o = — 2 + ^ + (T and po = — 4 + 2cr. 

Since the terms that grow stronger are those that involve a^, we shall restrict our study just to 
these. We shall try to obtain the power spectrum in a power-law form 

Vk{t) = C{t)k^''-^. (B.2) 

This will automatically lead us to a scale dependent spectral index rife, that it is found to be 



In 



rife = 3 - ^0 + 



In(fc) 



(B.3) 



- 16 - 



where C{t) = oqc 



-00 Hot 



depends on the first coefficient and the first power. We may write 

nk = na + Suk, 

no = 3 - /3o = 2/i ~ 3, 



In 



1 + E.=if^(fce"^°*)^°+' 



In(fc) 



(B.4) 
(B.5) 

(B.6) 



Since the values of k are related to super Hubble wavelengths : < k < 9ae^°^ (and assuming that 
(?cr <C 1), we see that < fce~^°* <C 1 and therefore it is pertinent a perturbative analysis in powers 
of 9a. In this case the dominant spectral index comes from uq — 2fi, and Sn^ are perturbative 



corrections. The integration of any of the power spectrums (B.l) or (B.2) provide us the amplitude 
for electric fields 



ekn 



dk 



Vk{t) 



„2Hot 



ao{9a) 



2+/3o 



ai{9a) 



3+Po 



1 4+^0 



2 + /?o 



3 + /3o 



4 + /3o 



(B.7) 



If we only stay with ng, this would mean we are cutting the previous expression just to the first 
term and only the coefficient Uq will appear. But if we keep to first order corrections in 6nk, we can 
see that the factor ai/ap appears in the correction. In general we shall obtain that to A^th-order 
correction, the first N coefficients will appear to each order respectively. 

In what follows we shall fix the spectral indices of the inflaton as /i = 3/2 + e, where e = — m^/2 



and TO^ is associated to the measured spectral scalar index tis ~ 0.96 pi 



^0.04. The 



spectral index of the vector fields is fixed so as to give a nearly scale invariant spectrum of the 
magnetic fields a = 5/2 + rj, with rj = —v^j^. For a similar spectrum to whole of the infiaton field 
it is expected to to be negative, but \v'^ \ <C 1. 

Studying just the terms that grow faster, as e^^°*, and considering sufficiently large scales 
9 <SC^ 1, we shall cut the power series to the first two terms. For the previous values of /i and <j there 
is no fli, since the power series begins after the aq term, in ai- 



Since (3/2 - = 2f to obtain ao only we need to find in ( IATLII ) 



00 



P2 i07r3/2 



and then 



ao 



27r225 



2/5 



(B. 



(B.9) 



Since both, e and r\ are respectively small departures from /i and cr, then we obtain that (e — r])^^ 3> 
1. We notice that 2/5(e — r;)^^ comes from the solution Xg, that considers only the contribution 
from the inhomogeneous solution of ^ coupled to the effective infiaton. This means that here the 
most relevant solution is , and only considering this solution one obtains 



ao = 



2ij4 



54 (e - ^f 



(B.IO) 



with a spectral index 710 = 3 and buk = 0. 
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Figure 1: The figure shows 57ik as a function of ^ = k/kn- Notice that when the horizon entry (i.e., after 
inflation when 9 = 9, = k-j^/kn, ~ 1), the value of Suk is close to Sn^, — —0.035. In all our calculations we 
have taken into account the value kn, = (re®" as the wavenumber related to the horizon wavelength when, 
after inflation, the horizon entry. 
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